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Abstract: We present the analysis and results of 12.5 hours of high-energy gamma-ray observations 
of the EGRET-detected pulsar PSR B 195 1+32 using the Solar Tower Atmospheric Cherenkov Effect 
Experiment (STACEE). STACEE is an atmospheric Cherenkov detector, in Albuquerque, New Mexico, 
that detects cosmic gamma rays using the shower-front-sampling technique. STACEE's sensitivity to 
astrophysical sources at energies around 100 GeV allows it to investigate emission from gamma- ray 
pulsars with expected pulsed emission cutoffs below 100 GeV. We discuss the observations and analysis 
of STACEE's PSR 1951+32 data, accumulated during the 2005 and 2006 observing seasons. 



Introduction 

Young energetic pulsars (rapidly rotating, highly 
magnetized neutron stars, that produce non- 
thermal photons) remain as yet the only Galac- 
tic objects unambiguously detected within the 
EGRET energy range. To date, more than 1500 
pulsars have been observed at radio energies [1]. 
About 70 of these are seen in X-rays [2] but only 6 
(Crab, Geminga, Vela, PSR B 195 1+32, PSR1706- 
44, and PSR B 1055-52) were detected by EGRET 
[3]. Of these, only one, PSR B1951+32, has been 
observed up to ^20 GeV without any apparent cut- 
off in its differential energy spectrum [4]. Accord- 
ingly, PSR 1951+32 is often considered the best 
pulsar candidate for pulsed-emission detection in 



the very-high-energy (VHE) regime; i.e. at ener- 
gies above ~50 GeV. 

As is the case for most pulsars, PSR B1951+32 
was first detected at radio energies. It was discov- 
ered with a 39.5 ms period in the radio synchrotron 
nebula CTB 80 [5]. From the radio observations 
it was deduced that the pulsar has a characteristic 
age of l.lxlO 5 yr with a surface magnetic field 
of 4.9 xlO 11 G, and a rotational energy loss rate 
of 3.7xl0 36 ergs s _1 . It has also been observed 
in X-rays. While the pulsar emits a single pulse at 
radio frequencies and in X-rays, the EGRET obser- 
vations display a double-pulsed profile with neither 
of the gamma-ray peaks coinciding with the radio 
peak. Interestingly, there is no evidence for any 
interpeak emission in the gamma-ray data [4]. 
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At TeV energies, only upper limits on the pulsed 
emission, and on the emission from the surround- 
ing synchrotron nebula, exist [6]. A recent report 
by the MAGIC collaboration [7] presents observa- 
tions above 75 GeV with no evidence for pulsed 
emission. 

Models that attempt to explain the non-thermal 
high-energy pulsed emission from pulsars gener- 
ally fall into one of two broad categories; the Po- 
lar Cap [8] or Outer Gap [9] models. While both 
models can explain the observed gamma-ray emis- 
sion at EGRET energies, they differ in their pre- 
dictions for detectable emission above 20 GeV. In 
both scenarios the pulsed emission is attributed to 
particle acceleration in the pulsar's magnetosphere, 
with spectral cutoffs predicted in the 20-100 GeV 
energy range. 

The Polar Cap model localizes the emission site 
to a region close to the magnetic poles of the neu- 
tron star, where the magnetic field is strong. The 
Outer Gap model, on the other hand, contends that 
gamma-ray production occurs far from the neutron 
star surface in a region of relatively weak magnetic 
field, in so called "outer gaps" near the null surface 
of the outer magnetosphere. 

The emission sites dictate the energy of the ex- 
pected spectral cutoff, insofar as the maximum en- 
ergy of the curvature-radiated photons escaping the 
magnetosphere is limited by pair-production in the 
pulsar's magnetic field. Since the magnetic field is 
stronger near the polar cap, the Polar Cap model 
anticipates a lower-energy cutoff than the Outer 
Gap model. 

Any detection of TeV emission would clearly fa- 
vor the Outer Gap model and would thereby sig- 
nificantly contribute to our understanding of pul- 
sar emission processes. Thus far, no pulsed TeV 
emission has been detected from any pulsar, al- 
though ever-lower upper limits are constraining 
the emission models. As a gamma-ray observa- 
tory operating at the lower end of the VHE regime 
(around 100 GeV) the Solar Tower Atmospheric 
Cherenkov Effect Experiment (STACEE) is suited 
to pulsar observations. STACEE observations of 
PSR B 195 1+32, undertaken in 2005 and 2006, are 
reported here. 



STACEE 

STACEE is a shower-front-sampling atmospheric 
Cherenkov telescope that uses the facilities of the 
National Solar Thermal Test Facility (NSTTF) in 
Albuquerque, New Mexico [10]. The NSTTF is a 
solar energy research facility incorporating a cen- 
tral receiver tower and an array of heliostats (so- 
lar mirrors). STACEE uses secondary mirrors, in 
the central receiver tower, to focus Cherenkov light 
from air showers, that is reflected by the heliostats, 
onto cameras having a total of 64 photomultiplier 
tubes (PMTs). A one-to-one mapping between he- 
liostats and PMTs allows the Cherenkov shower- 
front to be sampled independently at 64 different 
locations, spread over an area of ^2xl0 4 m 2 , on 
the heliostat field. The large reflecting surface pro- 
vided by the 64 heliostats, each with an area of 
37 m 2 , allows STACEE to operate with an energy 
threshold around 100 GeV. 

STACEE uses a custom-built 2-level trigger sys- 
tem [10] to select Cherenkov events from amongst 
the background of night-sky light fluctuations. In 
the event of a Cherenkov trigger, amplified and 
AC-coupled signals from the PMTs are recorded, 
together with a GPS timestamp (1/xs resolution), 
using 8-bit Flash Analog to Digital Convertors 
(FADCs), one per PMT The FADCs provide im- 
portant temporal and intensity information, at a 
sampling rate of 1 GS/s, which is fully utilized in 
the offline data analysis procedure. 

Observations of PSR B1951+32 

STACEE observations of PSR B 195 1+32 were un- 
dertaken during clear moonless nights in the pe- 
riods June-July 2005, September 2005 and June 
2006. The declination of PSR B1951+32 is favor- 
able for high-elevation observations by STACEE 
and most data were recorded at an elevation of 82.5 
degrees. 

Observation runs of 28 minute duration were con- 
ducted in an on-source-only mode. This is differ- 
ent from STACEE's normal "ON/OFF" observing 
mode in which observations of both the source re- 
gion (ON) and a control region of sky (OFF) are 
recorded. For PSR B 195 1+32, on-source-only ob- 
servations were considered appropriate, given the 
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lack of interpeak gamma-ray emission observed 
by EGRET [4] and the desire to maximize the 
amount of data recorded. Should the source emit 
gamma rays above 100 GeV within the sensi- 
tivity of STACEE, they would be detectable as 
an excess number of events within the EGRET 
pulsed region (phases ranges 0.12-0.22 and 0.48- 
0.74), compared to the number of events outside 
the pulsed region, assuming that the EGRET pulse 
profile is maintained at TeV energies. 

A total of 15.1 hours of PSR B 195 1+32 data were 
recorded by STACEE during 2005 and 2006. 

Data Quality Selection 

Before STACEE data are analyzed for the pres- 
ence of a gamma-ray signal, data quality selection 
is performed to remove sections of data flagged 
as unusable due to hardware malfunctions, unsta- 
ble atmospheric conditions or transient terrestrial 
light contamination. Heliostat and FADC mal- 
functions are logged each night by the heliostat 
and FADC control software respectively, and later 
merged with the Cherenkov data offline. Hence, 
data contaminated by hardware malfunctions are 
easily identified and removed. 

Data contaminated by unstable atmospheric con- 
ditions or transient background light are usually 
identified by cross-checking the level- 1 trigger 
rates (the rates at which clusters of 8 heliostats reg- 
ister coincident threshold-crossing light pulses) be- 
tween ON and OFF data. Under stable conditions, 
the ON and OFF level- 1 trigger rates should be- 
have similarly. For the PSR B 195 1+32 observa- 
tions reported here, a modified data quality check- 
ing procedure was necessary since OFF data were 
not available for comparison. The following pro- 
cedure was adopted. 

For each data run recorded, the level- 1 trigger 
rates were binned in 10-second time slices and 
lightcurves of the rates as a function of time 
produced. Second-order polynomial fits to the 
lightcurves were used to characterize the elevation- 
dependence of the trigger rates over the course of 
the data run. The residual difference between the 
true rate and the fitted rate in each time slice was 
then examined and slices with a residual more than 
3 times the RMS residual were flagged as bad. 



Time slices for which three or more clusters were 
labeled as bad were cut from the data set, together 
with 25-second adjacent time sections. In this way, 
1 minute periods of data with anomalous level- 1 
trigger rates — indicative of unstable weather or lo- 
cal light contamination — were excluded from the 
final data analysis. 

Data quality selection removed approximately 2.6 
hours of PSR B 195 1+32 data, such that the final 
data set used in the gamma-ray and temporal anal- 
yses was 12.5 hours of clean data. 

Data Analysis 

Data analysis for STACEE is a two-stage process. 
In stage 1, fixed and dynamic time corrections 
are applied to the pulse times registered for each 
PMT Fixed delays are necessary to account for 
the measured photon time-of-flight delays from the 
individual heliostats to the PMTs and for the sig- 
nal propagation delays between the PMTs and the 
central trigger system. Dynamic delays account 
for the different and changing arrival times of the 
Cherenkov shower-front at the heliostats. 

In stage 2, the time-corrected FADC data are used 
to determine the shower core — the point at which 
the shower's axis intersects the Earth's surface. A 
template-fitting approached is used. PMT-charge 
templates are compiled using showers simulated 
over a large range of zenith angles, azimuth an- 
gles, and core locations. By finding the template 
that best matches a particular event, an approxi- 
mate core location for that event is obtained. Core 
resolution using this method is about 15 m, accord- 
ing to simulations. 

The estimated shower core location is used as 
the starting point of the "grid ratio" analysis — 
the analysis routine used by STACEE for 
gamma/hadron discrimination [11, 12]. The grid 
ratio essentially provides a measure of the spheric- 
ity of the shower-front. Since gamma-ray show- 
ers with energies in the STACEE operating range 
typically have spherical shower fronts, whereas 
hadron showers are irregular and non-spherical, a 
cut on the sphericity of the shower-front allows for 
a significant reduction of the large hadronic back- 
ground. 
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Figure 1: Phaseogram of candidate STACEE 
gamma-ray events for PSR B1951+32 folded us- 
ing the pulsar's radio ephemeris. The shaded re- 
gions show the positions of the main pulse (phase 
0. 12-0.22), and the interpulse (phase 0.48-0.74) as 
seen by EGRET [4] 



cess events within the pulsed phase profile seen 
by EGRET at lower energies [4]. That is, emis- 
sion is assumed to occur in the phase range of the 
main pulse, phase 0.12-0.22, and the interpulse, 
phase 0.48-0.74. A differential flux upper limit of 
4.52 x 10~ 6 MeV cm~ 2 s _1 was determined at the 
energy threshold of 117 GeV, by extrapolation of 
the EGRET spectrum to STACEE energies. 
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Temporal Analysis 

The arrival times of the candidate gamma-ray 
events, as registered by the STACEE GPS clock, 
were transformed to the solar system barycen- 
ter using the JPL DE200 Planetary and Lunar 
Ephemeris. Two independent barycentering algo- 
rithms were employed with excellent agreement. 
The pulsar phase of each event at the solar system 
barycenter, with respect to the PSR B 195 1+32 ra- 
dio ephemeris, was calculated and a phaseogram 
produced. 

To test the barycentering algorithms, an optical sig- 
nal of the Crab pulsar, as seen using three PMTs 
set up at STACEE in a special arrangement, was 
analyzed. The resulting optical lightcurve, showed 
clear evidence of a pulsed signal phase-aligned 
with the Crab pulsar radio pulsations [13]. This 
demonstrated the validity of the barycentering soft- 
ware and also STACEE's timing electronics. 

Results and Conclusions 

No evidence for pulsed gamma-ray emission above 
117 GeV was found in the selected gamma-ray 
data set used in this work (Figure 1). To calcu- 
late flux upper limits, we used the method of He- 
lene [14] to estimate 99.0% upper limits for ex- 
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